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I. INTRODUCTION 
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A general analytic procedure is developed to deal with the Newtonian limit of f(R) gravity. A 
discussion comparing the Newtonian and the post-Newtonian limit of these models is proposed in 
order to point out the differences between the two approaches. We calculate the post-Newtonian 
parameters of such theories without any redefinition of the degrees of freedom, in particular, without 
adopting some scalar fields and without any change from Jordan to Einstein frame. Considering the 
Taylor expansion of a generic f(R) theory, it is possible to obtain general solutions in term of the 
metric coefficients up to the third order of approximation. In particular, the solution relative to the 
gu component gives a gravitational potential always corrected with respect to the Newtonian one 
of the linear theory f(R) = R. Furthermore, we show that the Birkhoff theorem is not a general 
result for /(ii)-gravity since time-dependent evolution for spherically symmetric solutions can be 
achieved depending on the order of perturbations. Finally, we discuss the post-Minkowskian limit 
and the emergence of massive gravitational wave solutions. 

PACS numbers: 04.25.-g; 04.25. Nx; 04.40.Nr 

cr 

In recent years, the effort to give a physical explanation to the today observed cosmic acceleration [l|, 0, Q has 
, attracted a good amount of interest in /(i?)-gravity, considered as a viable mechanism to explain the cosmic accel- 
eration by extending the geometric sector of field equations 0, S @ • There are several physical and mathematical 
motivations to enlarge General Relativity (GR) by these theories. For comprehensive review, see B, M, B- 

Specifically, cosmological models coming from /(i?)-gravity were firstly introduced by Starobinsky [lfj in the early 
80'ies to build up a feasible inflationary model where geometric degrees of freedom had the role of the scalar field 
ruling the inflation and the structure formation. 

On the other side, dealing with such extended gravity models at shorter astrophysical scales (Galaxy and Solar- 
System), one faces the emergence of corrected gravitational potentials with respect to the Newton one coming out 
from GR. This result is well known since a long time [ll[ , and recently it has been pursued to carry out the possibility 
to explain the flatness of spiral galaxies rotation curves without the addition of huge amount of Dark Matter. In 
particular, the rotation curves of a wide sample of low-surface-brightness spiral galaxies have been successfully fitted 
by these corrected potentials [l2j and reliable results are also expected for other galaxy- types [l3l |. 

Other issues as, for example, the observed Pioneer anomaly problem [14| can be framed into the same approach 
(TH and then, apart the cosmological dynamics, a systematic analysis of such theories urges at short scale and in the 
low energy limit. 

In this paper, we are going to discuss, without specifying the form of the theory, the Newtonian limit of /(i?)-gravity 
pointing out the differences and the relations with respect the post-Newtonian and the post-Minkowskian limits. In 
literature, there are several definitions and several claims in this direction but clear statements and discussion on 
these approaches urge in order to find out definite results to be tested by experiments [20l |. 

The discussion about the short scale behavior of higher order gravity has been quite vivacious in the last years 
since GR shows is best predictions just at the Solar System level. As matter of fact, measurements coming from 
weak field limit tests like the bending of light, the perihelion shift of planets, frame dragging experiments represent 
inescapable tests for whatever theory of gravity. Actually, in our opinion, there are sufficient theoretical predictions 
to state that higher order theories of gravity can be compatible with Newtonian and post-Newtonian prescriptions. 



* e-mail address: capozziello@na.infn.it 
t e - mail address: stabile@sa.infn.it 
+ e-mail address: antrois@gmail.com 



2 



In other papers pjj, we have shown that this result can be achieved by means of the analogy of f(R)- models with 
scalar - tensor gravity. 

Nevertheless, up to now, the discussion on the weak field limit of f(R) - theories is far to be definitive and there are 
several papers claiming for opposite results (l8l . Il9| , or stating that no progress has been reached in the last forty due 
to the several common misconceptions in the various theories of gravity [20(. 

In particular, people approached the weak limit issue following different schemes and developing different parame- 
terizations which, in some cases, turn out to be not necessarily correct. 

The purpose is to take part to the debate, building up a rigorous formalism which deals with the formal definition 
of weak field and small velocities limit applied to fourth-order gravity. In a series of papers, our aim is to pursue a 
systematic discussion involving: i) the Newtonian limit of /(i?)-gravity (the present paper), ii) spherically symmetric 
solutions vs. the weak filed limit in f(R) -gravity 12211; and, finally, Hi) general fourth-order theories where also 
invariants as R^R^" or i? aj g At „i? a/3 ' 1I/ are considered, [2l| . 

Our analysis is based on the metric approach, developed in the Jordan frame, assuming that the observations are 
performed in it, without resorting to any conformal transformation as done in several cases [Hjj . This point of view 
is adopted in order to avoid dangerous variable changes which could compromise the correct physical interpretation 
of the results. 

We will show that the corrections induced on the gravitational potentials can be suitable to explain relevant 
astrophysical behaviors or can be related with some relevant physical issues. 

As a preliminary analysis, we will concentrate on the vacuum case with the aim to build up a further rigorous 
formalism for the Newtonian and post-Newtonian limit of f(R) theories in presence of matter. As we will see, it is 
possible to deduce an effective estimation of the post- Newtonian parameter 7 by considering the second order solutions 
of the metric coefficient in the vacuum case. For the sake of completeness we will treat the problem also by imposing 
the harmonic gauge on the field equations. 

The paper is organized as follows: in Sec. II, the general formalism concerning the spherically symmetric background in 
fourth order gravity is introduced; Sec. Ill is devoted to a discussion of the post-Newtonian approximation considering 
the differences with respect GR: in this theory not all order of perturbations can be consistently achieved if conservation 
laws are taken into account, in /(i?)-gravity this shortcoming can be, in principle, avoided. In Sec. IV, the analytic 
approach to the weak field in /(i?)-gravity is developed. In particular, we achieve the gravitational potential (related to 
the gtt-component of the metric) which is always corrected with respect to the Newtonian one of the linear f(R) = R 
theory. Besides, we show that the Birkhoff theorem is not a general result for /(i?)-gravity since time-dependent 
evolution for spherically symmetric solutions can be achieved depending on the order of perturbations. In Sec. V, the 
post-Minkowskian limit is discussed considering also the possibility to obtain gravitational waves solutions. Sec. VI 
is devoted to the discussion and conclusions. 

II. f(R) - GRAVITY IN SPHERICALLY SYMMETRIC SPACETIME 

The action for f(R)- gravity reads: 



A 



f(R) + XC T , 



(1) 



where f(R) is an analytic function of Ricci scalar, X = 16 ^ G is the coupling constant and C m describes the ordinary 
matter Lagrangian. Such an action is the straightforward generalization of the Hilbert-Einstein action of GR where 
f(R) = R is assumed. 

By varying (fT]) with respect to the metric, one obtains the fourth-order field equations : 

1 X 
f'Rpu - g/SV ~ f'fiu + 9^0 f = —T MV , (2) 



with T uv = ^^OAn) d f , = 4fW The trace is 



3Df' + f'R-2f= yT, (3) 

and such an expression can be read as a Klein-Gordon equation, where the effective field is /', if f(R) is non-linear 
in R [T(J. 
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As said, we are interested in investigating the Newtonian and the post-Newtonian limit of /[i?)-gravity in a 
spherically symmetric background. Solutions can be obtained considering the metric (see also [23l. |24|) : 

ds 2 = g aT dx a dx T = A(x°, r)dx° 2 - B(x°, r)dr 2 - r 2 dfl (4) 

where x° — ct; A and B are generic functions depending on time and coordinate radius; dfl is the angular element. 
The field equations (J2j) turn out to be 



R-fiv f Rfj,L> 2 f dfJ-v ^■fj,v 2 

H = g^H aT = f'R -2f + H=§T 



where 



(5) 



g w i0 + 5 UU In j—g A R + g r \ r + g rr In R r + 



+g 00 RflO + 9 rT R rr 

H = g° T H aT = if" 



-f" 



R„R,-9^[g 00 Ro 2 + g rr Rr 2 



f° + .g 00 In ^=5 o R fi + 9 rr , r + 9 rr ^ V=9. r ) R ,r + 9°°R,oo + g rr R,, 



(6) 



-if" 



g 00 R/ +g rr R/ 



are the higher than second order terms of the theory. We are adopting the convention R^ v = R p ^ pv for the Ricci 
tensor and R a p tlL , — rj| — for the Riemann tensor. Connections are Levi-Civita : 



r a/3 = ^5 mp (5«p,/3 + gp P , a - g a p,p) ■ 



(7) 



III. GENERAL REMARKS ON THE NEWTONIAN AND THE POST-NEWTONIAN 

APPROXIMATION 



At this point, it is worth discussing some general issues on the Newtonian and post-Newtonian limits. Basically 
there are some general features one has to take into account when approaching these limits, whatever the underlying 
theory of gravitation is. 

If one consider a system of gravitationally interacting particles of mass M, the kinetic energy ^Mv 2 will be, roughly, 
of the same order of magnitude as the typical potential energy U = GM 2 /f, with M, f, and v the typical average 
values of masses, separations, and velocities of these particles. As a consequence: 



(for instance, a test particle in a circular orbit of radius r about a central mass M will have velocity v given in 
Newtonian mechanics by the exact formula v 2 — GM/r.) 

The post- Newtonian approximation can be described as a method for obtaining the motion of the system to an higher 
than the first order (approximation which coincides with the Newtonian mechanics) with respect to the quantities 
GM I? and v 2 assumed small with respect to the squared light speed c 2 . This approximation is sometimes referred to 
as an expansion in inverse powers of the light speed. 

The typical values of the Newtonian gravitational potential U are nowhere larger than 10~ 5 in the Solar System (in 
geometrized units, U/c 2 is dimensionless) . On the other hand, planetary velocities satisfy the condition v 2 < U, whilc§ 



8 We consider here on the velocity v in units of the light speed c. 
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the matter pressure p experienced inside the Sun and the planets is generally smaller than the matter gravitational 
energy density pU, in other words * p/p < U. Furthermore one must consider that even other forms of energy in the 
Solar System (compressional energy, radiation, thermal energy, etc.) have small intensities and the specific energy 
density IT (the ratio of the energy density to the rest-mass density) is related to U by II < U (II is ~ 10~ 5 in the Sun 
and ~ 10~ 9 in the Earth [25j] ) . As matter of fact, one can consider that these quantities, as function of the velocity, 
give second order contributions : 



u ~v 2 ~ p/p ~n~o(2). 



(9) 



Therefore, the velocity v gives 0(1) terms in the velocity expansions, U 2 is of order 0(4), Uv of 0(3), UT1 is of 0(4), 
and so on. Considering these approximations, one has 



and 



d 



v V. 



(10) 



\d/dx°\ 



0(1). 



(11) 



Now, particles move along geodesies : 



which can be written in details as 



cPx % 
dx^ 2 



2P n 



dx r ' 



dx° 



f^l+rn dx ° dxT - n 
ds 2 ar ds ds 



dx m dx n 
'' dx° dx° 



2r" 



dx r ' 



dx° 



dx m dx n 
mn dx° dx° 



dx l 
dx^ 



(12) 



(13) 



In the Newtonian approximation, that is vanishingly small velocities and only first-order terms in the difference 
between g^ v and the Minkowski metric 77^, one obtains that the particle motion equations reduce to the standard 
result : 



d 2 x l 
dx° 2 



00 



1 dgpo 

2 dx l 



(14) 



d 2 x i 
dx^- 



to the order GM/r 2 , that 



The quantity 1 — goo is of order GM/r, so that the Newtonian approximation gives 
is, to the order v 2 /r. As a consequence if we would like to search for the post-Newtonian approximation, we need to 



compute 



dx°< 



to the order w 4 /r. Due to the Equivalence Principle and the differentiability of spacetime manifold, 



we expect that it should be possible to find out a coordinate system in which the metric tensor is nearly equal to the 



Minkowski one 77^, the correction being expandable in powers of GM/r 
metric developed as follows : 



In other words one has to consider the 



' ff0 o(z°,x) ~ 1 +^ 2 V,x) +g$(x°,x) + 0(6) 
ff0l ( 2 ; ,x)~^ ) (x ,x) + O(5) 
^•( a; ,x)~-5 ii +4 2) ( a; ,x) + O(4) 



(15) 



t Typical values of p/p are ~ 10 5 in the Sun and ~ 10 10 in the Earth [25 
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where 6{j is the Kronecker delta, and for the controvariant form of g^ v , one has 

' g 00 (x°,x) ~ 1 + gW°°(x°,x) +gW° (x°,x) + 0(6) 
< <? 0l (x°,x) ~ 5 ( 3 )° l (x°,x) + 0(5) 
^ g lj (x°,x) ~ + . 9 ( 2 )^'(x°,x) + 0(4) . 



(16) 



In evaluating we must take into account that the scale of distance and time, in our systems, are respectively set 
by f and f/v, thus the space and time derivatives should be regarded as being of order 



9 1 d v 

dx l f ' da; f 
Using the above approximations ([15]), P^|) , we have, from the definition ([7]), 



(17) 



r(3)° - I„( 2 )'° 

1 00 — 2 WOO 



r(3) 1 -i( (3)' 4 
1 oj - 2 y oj 



' oj 9 j,o 



r( 4 ) 1 = l( a (*)L + ( 2 ) im o (2) -2o< 3 ) 1 

1 00 — 2 I # 00 ' i/00,m z i/ 0,0 



The Ricci tensor component are 



1 oo — 2-yoo 



r( 3 ) ■ = if ( 3 ' u ■ + ( 3 ' u ■ - 



r (4) o 4 -^ff (4) l + 5 (2)00 ^ 
r (2) o» = ^ (2) o. 



(18) 



K-( 4 ) - IA„< 4 ) _ I„(2) mn „(2) _ I„(2)" m „(2) , 1 (2) m _ 1 (2)°>"» (2) _ 1 „(2) m >" ( 2) _ (3) m 
n 00 — 2^500 29 ,m 9oo, n 2^ 9o0,mn ^ 2# m,00 4# SoO,™ 4^ m #00, ri 5 0,', 



(2) _ 
00 — 




(4) 

00 " 




(3) _ 
Oi — 




(2) 
ij 





m.Oi 



(19) 



and assuming the harmonic gauge g^T^ = (see the Appendix for details), one can rewrite these last expressions 
as 



f i? (2) - 


|A^ 


-^oo 


|A<$ - ^ 






R {2) = 





) -I (2)° _I|V7 ( 2 )|2 
Kmn 29 0,00 2 I Vi) J/QO I 



with A and Vj respectively, the Laplacian and the gradient in flat space. The Ricci scalar reads 



(20) 



i?( 2 ) = R(V° - R^Z = IA 5 ( 2 )° - \£>gW 



R(4) _ p(4)° , „(2) 00 7 ? (2) , „(2) m "p(2) _ 1 a (4)° 1„(2)°'° 1 _(2) 
H - K + .9 ^00 + ' -Kmn - 2 Aff V ; _ 2~9 0,0 ~ 2"S 



5 (2) o 



Ao (2) 

{—^■ymn 



ilV ff (2) ol 2 + ^ (2)0 °A^ 
(21) 
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The inverse of the metric tensor is defined by means of the equation 



= dp 



(22) 



with #g the Kronecker delta. The relations among the higher than first order terms turn out to be 



,(4)00/ 



(x , x) = g$ (x , x) - g^l> (x , x) 



I 5 (2) ^o,x) 



(23) 



7|?(^o,x) 



Finally the Lagrangian of a particle in presence of a gravitational field can be expressed as proportional to the invariant 
distance ds 1 / 2 , thus we have: 



dx p dx a 



1/2 



V2 =U + 9$ + 9® + 2g^v m v 2 + gWv m v n ) ^ , (24) 



= I # + ZgOmV™ + g,nnV m V 

which, to the 0(2) order, reduces to the classic Newtonian Lagrangian of a test particle LNew 



(2) 

•Soo 



1/2 



where v — ^Tr^rr- As matter of fact, post-Newtonian physics has to involve higher than 0(4) order terms in the 
Lagrangian. 

An important remark concerns the odd-order perturbation terms 0(1) or 0(3). Since, these terms contain odd 
powers of velocity v or of time derivatives, they are related to the energy dissipation or absorption by the system. 
Nevertheless, the mass-energy conservation prevents the energy and mass losses and, as a consequence, prevents, in the 
Newtonian limit, terms of 0(1) and 0(3) orders in the Lagrangian. If one takes into account contributions higher than 
0(4) order, different theories give different predictions. GR, for example, due to the conservation of post-Newtonian 
energy, forbids terms of 0(5) order; on the other hand, terms of 0(7) order can appear and are related to the energy 
lost by means of the gravitational radiation. 



IV. THE NEWTONIAN LIMIT OF f(R) GRAVITY IN SPHERICALLY SYMMETRIC BACKGROUND 

VS. POST-NEWTONIAN LIMIT 

Exploiting the formalism of post-Newtonian approximation described in the previous section, we can develop a 
systematic analysis in the limit of weak field and small velocities for /(i?)-gravity. We are going to assume, as 
background, a spherically symmetric spacetime and we are going to investigate the vacuum case. Considering the 
metric ((3]), assuming, unless not specified, c — 1 and then x° = ct — > t, we have, for a given g^ : 



9tt(t,r) : 


= A(t,r)~ 


9rr(t,r) 

< 


= -B(t,r) 


9ee(t,r) 




> 94>4>^ r ) 


9 • 2 

= — r sm 



,(2), 



.(4)/ 



while the approximations for g pv are 



.(2)/ 



(25) 



7" = A(^r)- 1 ~ 



1 - 



(2) 

9tt 



[9u ] ~9u 



(4)l 



-l-o (2) 

l y rr 



(26) 
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The determinant reads 



g ~ r 4 sin 2 + [g$ - 3t (2) ] + [St ( M } - fl& } ]} ■ 
As a consequence, the Christoffel's symbols are 



(27) 



* (2) 

r(3)* _ »k± 

L tt — 2 



r(3)' 



r(3)* = 

rr 2 



= sin 2 0r£ 



1 tt - * -1 tt — 2 2 



(2) 



(4) _ (2) (2) 



T (2V r , r (4)^ _ 9% r alVaJil 



r(o) r + r( 2 ) r +rW r -- r -™ (2) -™ (2) 



Let us even display the Ricci's tensor components 



where 



Rtt 

Rtr 
Rrr 



d(2) , n(4) 



p(3) 



Ree — R, 



(2) 



i?^ ~ sin 6R 



2/) R (2) 



(28) 



(29) 



(2) _ -rg^> rr +2g t ( t 2) r 



R 



R 



(3) 



i? 



(4) _ -r^^'+4 g {tU^L 2 U^ r +2g^H2 g L 2 U^L 2 ^]+2r- 3 i^ r +2r g 



4r 



(2) _ ^3j t ; rr +2g 



— — 
1 Ltv — 



< 2 » 4-9„( 2 ) 



2r 



(2) _ 2 fl ff+r[ g < f a l+ g ff,,] 



(30) 



and the Ricci scalar expression in the post-Newtonian approximation 



R ~ i? (2) + i? (4) 



with 



(31) 



i?( 2 ) 



2 g < 2 ;+W2 g < 2 ,t+2 g ( 2 r ) r +^ t (2 ,t.] 



(32) 



R (4) = ^r+^V[2g[i] r +i g ^ r+ r g [i] r ^ 



In order to derive the post-Newtonian approximation for a generic function f(R), one should specify the f(R)- 
Lagrangian into the field equations ([5]). This is a crucial point because once a certain Lagrangian is chosen, one 
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will obtain a particular post-Newtonian approximation referred to such a choice. This means to lose any general 
prescription and to obtain corrections to the Newtonian potential which refer "univocally" to the considered f(R) 
function. Alternatively, one can restrict to analytic f(R) functions expandable with respect to a certain value R = Rq. 
In general, such theories are physically interesting and allow to recover the GR results and the correct boundary and 
asymptotic conditions. Then we assume 



f(R) = E L± F Z CR - RoT ^fo + fiR + hR 2 + f 3 R 3 + ... • (33) 

n 

On the other hand, it is possible to obtain the post-Newtonian approximation of /(i?)-gravity considering such an 
expansion into the field equations © and expanding the system up to the orders O(0), 0(2) e 0(4). This approach 
provides general results and specific (analytic) Lagrangians are selected by the coefficients fi in (|33|) . 

Let us now substitute the series (13"3"|) into the field Eqs. ((SJ. Developing the equations up to 0(0), 0(2) and 0(4) 
orders in the case of vanishing matter, i.e. T^ v = 0, we have 



rr(0) 


= o, 


H<f» 


= 


rr(2) 


= 0, 




= 


rr(3) 


= o, 


ff(3) 


= 


rr(4) 


= 0, 


H w 


= 



(34) 



and, in particular, from the O(0) order approximation, one obtains 



fo = . 



(35) 



This result suggests a first consideration. If the Lagrangian is developable around a vanishing value of the Ricci scalar 
(Rq = 0) the relation (|35|) will imply that the cosmological constant contribution has to be zero in vacuum whatever 
is the /(-R)-gravity theory. This result appears quite obvious but sometime it is not considered in literature. 
If we now consider the 0(2)- order approximation, the equations system (|34p . in the vacuum case, results to be 



' hrR^ 


- 2/isgl + 


8/ 2 i?, ( r } 


f (2) 
Jl r 9tt,rr 


+ 4/ 2 ri?( 2 ) = 


ArRW 


Of „ (2) 4 
lyrr,r \ 


- 8/ 2 i?^ 2) 


f (2) 


= 


Of 


- r[firR& 


f n {2) 
- Jl9tt,r 


f n {2) 4- 
J iyrr,r ~r 


4/ 2 i?^ 2) + Af 2 rR$] = 


hrR^ 


+ 6/ 2 [2i?, ( r 2) 


+ rR.fr] 


= 





(36) 



<lg?) + r[2g^ r 



rR^+2g r % + rg^ rr }=0 



(2) 



The trace equation (the fourth line in the (|36p ). in particular, provides a differential equation with respect to the 
Ricci scalar which allows to solve the system (|3f)|) at 0(2)- order: 



(2) X 

9tt =o - 
J2) _ ftfflfr 



Si(t)e~ 



S 2 (t)e r 



3£r 1 



6(-?) 3/2 '- 



R (2) _ Si{t)e 



3£r 



6{ 2 r 



(37) 



2f r 



fl 

where £ = — — and /i and / 2 are the expansion coefficients obtained by Taylor developing the analytic f(R) La- 

6/2 

grangian. Let us notice that the integration constant So is correctly dimensionless, while the two arbitrary functions 
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of time 6x(t) and 5 2 (t) have respectively the dimensions of lenght^ 1 and lenght~ 2 ; £ has the dimension lenght~ 2 . 
The functions 5i(t) (i = 1,2) are completely arbitrary since the differential equation system (|36[) contains only spatial 
derivatives. Besides, the integration constant 5q can be set to zero, as in the theory of the potential, since it represents 
an unessential additive quantity. 

With these results in mind, the gravitational potential of a generic analytic f(R) can be obtained. In fact, the first 
of (I37|) gives the second order solution in term of the metric expansion (see the definition ([25]) ). but, as said above, this 
term coincides with the gravitational potential at the Newtonian order. In other words, we have gu = 1 + 20 gral) = 

(2) 

1 + g tt and then the gravitational potential of a fourth order gravity theory, analytic in the Ricci scalar R, is 



b FOG 
grav 



K ie - r ^ K 2 e r ^ 



3£r 



6(-0 



3/2 



(38) 



with K x = di(t) and K 2 = 6 2 (t). 

As first remark, one has to notice that the structure of the potential, for a given f(R) theory, is determined by 
the parameter £, which depends on the first and the second derivative of the f{R) function, once developed around a 
particular point Rq. 

Furthermore, one has to consider that the potential (|38[) holds in the case of non- vanishing f 2 since we manipulated 
the equations in (|36|) dividing by such a quantity. As matter of fact, the GR Newtonian limit cannot be achieved 
directly from the solution (|38|) but from the field equations (|36j) once the appropriate expressions in terms of the 
constants /,; are derived. 

The solution ([38]) has to be discussed in relation to the sign of the term under the square root in the exponents. 
The first possibility is that the sign is positive, which means that /i and f 2 have opposite signature. In this case, 
the solutions (|3"Tj) and can be rewritten introducing the scale parameter I = |£| -1 / 2 . In particular, considering 
Sr, = 0, the two Si(t) functions as constants, k\ = (5i(t)/3)l and k 2 (t) — (S 2 (t)/6)l 2 and by introducing a radial 
coordinate f in units of I, we have: 



9tt — °0 H a WT H r — 777 - K i— + K 2 — 



3 r/l 1 6 r/l 
(2) _ ( r /l+l)e- r /' , S 2 (t)l 2 (r/l-l)e^ 1 



g rr 



r/l 



R(2) _ <?i(«) e-"/' , 3 2 (t) e*-/' _ 3 
^ — I r/l "f" 2 r/J ~ F 



r/i 



-fcj 



(f+l)e 



+ /c 2 



(r-l)e 



(39) 



by which we can recast the gravitational potential as 



b FOG 
grav 



k 2 e r 



(40) 



which is analogous to the result in [ll(, derived for the theory R + aR 2 + (iR ilv R til ' and coherent^ with the results 
in Ref . [26l| . obtained for higher order Lagrangians as f(R,DR) = R + J2k=o akRO k R. In this last case, it was 
demonstrated that the number of Yukawa corrections to the gravitational potential was strictly related to the order 
of the theory. However, as discussed in [2l|, it is straightforward to show that the usual form Newton + Yukawa can 
be easily achieved by Eq. ([40|) through a coordinate change. 

From (|3T|) and (f3"9"]) . one can notice that the Newtonian limit of any analytic /(i?)-theory is related only to the first 
and second term of the Taylor expansion of the given theory. 

In other words, the gravitational potential is always characterized by the two Yukawa corrections and only the first 
two terms of the Taylor expansion of a generical f(R) Lagrangian turn out to be relevant. This is indeed a general 
result. 

The diverging contribution, arising from the exponential growing mode, has to be carefully analyzed and, in 
particular, the physical relevance of this term must be evaluated in relation to the length-scale (— £) _1 / 2 . For very 
large r, (i.e. r » (— _1 ^ 2 )i the weak field approximation turns out to be unphysical and the (|37|) does not hold 



^ Let us remember that in the case of homogeneous and isotropic spacetime, higher order curvature invariants as R I1V R ,1V and R a 
reduce to R 2 . 



,R a Pv- v 
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anymore. As matter of fact, one can obtain a modified gravitational potential which can work as a standard Newtonian 
one, in the opportune limit, and provide interesting behaviors at larger scales, even in presence of the growing mode, 
once the constants in the ([38]) have been opportunely adjusted. Such a potential, once the growing exponential term 
is settled to zero, reproduce the Yukawa-like gravitational potential, phenomenologically introduced by Sanders (27j 
to explain the flat rotation curves of spiral galaxies without dark matter. 

Besides, Yukawa-like corrections to the gravitational potential have been suggested in several approaches. For 
example, an interesting proposal is a model describing the gravitational interaction between dark matter and baryons. 
This points out that the interaction suppressed on small subgalactic scales can be described by means of a Yukawa 
contribution added to the standard Newtonian potential. Such a behavior is effectively suggested by the observations 
of the inner rotation curves of low-mass galaxies and provides a natural scenario in which to interpret the cuspy profile 
of dark matter halos observed in N-body simulations [28[ . 

It is important to stress that the result we have obtained here is coherent with other calculations. In fact, since 
the Taylor expansion of an exponential potential is a power law series, it is not surprising to obtain a power law 
correction to the Newtonian potential [12| when a less rigorous approach is considered in order to calculate the weak 
field limit of a generic /(i?)-theory. In particular, perturbative calculations will provide effective potentials which can 
be recovered by means of an appropriate approximation from the general case (|40p . 

Let us now consider now the opposite case in which the sign of £ is negative and, as a consequence, the two Yukawa 
corrections in (|39p are complex numbers. 

Since the form of g u , the gravitational potential (HO)) turns out to be: 



b FOG 
grav 



(41) 



which can be recast as 



b FOG 
grav 



— \(k\ H~ £2) cos r + z(&2 — fci) sin r] . 



(42) 



Such a gravitational potential, which could be discarded as a non-physically relevant, has the property to satisfy the 
Helmholtz equation, V 2 </> + k 2 (f> = AirGp, where <f> is the gravitational potential and p is a real function acting both 
as matter and the antimatter density. As discussed in [29| , Re [4>g^] can be addressed as a classically modified 
Newtonian potential corrected by a Yukawa factor while Im [<j)g^] could have significant implications for quantum 
mechanics. In particular, this term can provide an astrophysical, and in our case even theoretically well founded, 
origin for the puzzling decay Ki, — > tt + tt~ whose phase is related to an imaginary potential in the kaon mass matrix. 
Of course, these considerations, at this level, are only speculative, nevertheless it could be worth taking them into 
account for further investigations. 

Let us now consider the system (|34|) up to the third order contributions. The first important issue is that, at this 
order, one has to consider even the off-diagonal equation 



fi9 r % + 2f 2 rRf r = , (43) 

(2) 

which relate the time derivative of the Ricci scalar to the time derivative of g„ ■ . From this relation, it is possible 
to draw a relevant consideration. One can deduce that, if the Ricci scalar depends on time so it is for the metric 
components and even the gravitational potential turns out to be influenced. This result agrees with the analysis 
provided in where a complete description of the weak field limit of fourth order gravity has been provided in 
term of the dynamical evolution of the Ricci scalar. In that paper, it was demonstrated that if one supposes a time 
independent Ricci scalar, static spherically symmetric solutions are allowed. Eq. (|43p confirms this result and provides 
the formal theoretical explanation of such a behavior. In particular, together with the (|3"5|) . it suggests that if one 
considers the problem at a lower level of approximation (i.e. the second order) the background spacetime metric 
can have static solutions according to the Birkhoff theorem; this is no more verified when the problem is faced with 
approximations of higher order. In other words, the debated issue to prove the validity of the Birkhoff theorem in 
the higher order theories of gravity, finds here its physical answer. In [22| and here, the validity of this theorem is 
demonstrated for f(R) theories only when the Ricci scalar is time independent or, in addition, when the Newtonian 
limit solutions are investigated up to the second order of approximation in term of a 11/c expansion of the metric 
coefficients. Therefore, the Birkhoff theorem does not represent a general feature in the case of fourth order gravity 
but, on the other hand, in the limit of small velocities and weak fields (which is enough to deal with the Solar System 
gravitational experiments), one can assume that the gravitational potential is effectively time independent according 
to (J3TJ) and ([38]). 
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The above results fix a fundamental difference between GR and fourth order gravity theories. While in GR a 
spherically symmetric solution represents a stationary and static configuration difficult to be related to a cosmological 
background evolution, this is no more true in the case of higher order gravity. In the latter case, a spherically 
symmetric background can have time-dependent evolution together with the radial dependence. In this sense, a 
relation between a spherical solution and the cosmological Hubble flow can be easily achieved. 

The subsequent step concerns the analysis of the system ([34]) up to the 0(4) order. Such an analysis provides 

the solutions, in term of , the right order for the post-Newtonian parameters. Unfortunately, at this order of 
approximation, the system turns out to be too much involuted and a general solution is not possible. 

From Eqs. ([3"4l . one can notice that the general solution is characterized only by the first three orders of the f(R) 
expansion. Such a result is in agreement with the f(R) reconstruction which can be induced by the post-Newtonian 
parameters adopting a scalar-tensor analogy (for details see [171 130|). 

However, although we cannot achieve a complete description, an approximate estimation of the post-Newtonian 
parameter 7 can be obtained recurring to the 0(2) evaluation of the metric coefficients in the vacuum case. 

It is important to notice that, since (|37|) suggests a modified gravitational potential (with respect to the standard 
Newtonian one) as a general solution of analytic f(R) gravity models, there is no reason to ask for a post-Newtonian 
description for these theories. In fact, as previously said, the post-Newtonian analysis presupposes to evaluate 
deviations from the Newtonian potential at a higher than second order approximation in term of the quantity v/c. 
Thus, if the gravitational potential deduced from a given f(R) theory of gravity is a general function of the radial 
coordinate, displaying a Newtonian behavior only in a certain regime (or in a given range of the radial coordinate), 
it could be meaningless to develop a general post-Newtonian formalism as in GR [25|, [3l[. Of course, by a proper 
expansion of the gravitational potential for small values of the radial coordinate, and only in this limit, one can develop 
an analogous of the post-Newtonian limit for these theories with respect to the Newtonian behavior and estimate the 
deviations from it. 

In order to have an effective estimation of the post-Newtonian parameter 7, we can proceed in the following way. 
Expanding g tt and g rr , obtained at the second order in (f39|) with respect to the dimensionless coordinate f, one has** 



( 2 ) (u u \ 1 kl + fc2 1 kl + fc2 ;j nrol 
9tt = ( fc 2 - fci) + — ; — + — - — r + 0[2\ , 



(44) 



,(2) - 



ki + k 2 ki + k 2 _ 



-J rr 



-r - 



0[2] 



where, clearly, ki + k 2 = GM and ki — k 2 in the standard case. When f — * (i.e. when the coordinate r « \J— £) 
the linear and the higher than first order terms are vanishingly small and only the first Newtonian term survives. 
Since the post-Newtonian parameter 7 is strictly related to the coefficients of the 1/r term into the expressions of gtt 
and g rr , actually one can obtain an effective estimation of this quantity confronting the coefficients of the Newtonian 
terms relative to both the expressions in (|44[) . Being 7 = 1 in GR, the difference between these two coefficients gives 
the effective deviation from the GR expectation value. 

It is easy to derive that a generic fourth order gravity theory provides a post-Newtonian parameter 7 which is 
consistent with the GR prescription (7 = 1) if ki = k 2 . Conversely, deviations from such a behavior can be 
accommodated by tuning the relation between the two integration constants ki and k 2 . This is equivalent to adjust 
the form of the f(R) theory in such a way to obtain the right GR limit, and then the Newtonian potential. This 
result agrees with the viewpoint that asks for the recovering of GR behavior from generic f(R) theories in the post- 
Newtonian limit [32l . l33l ]. This is particularly true when the f(R) Lagrangian behaves, in the weak field and small 
velocities regime, as the Hilbert-Einstein Lagrangian. 

On the other side, if deviations from these regime are observed, a f(R) Lagrangian, built up with a third order 
polynomial in the Ricci scalar, can be suitable to interpret such a behavior (see [301]'). 

Actually, the degeneracy regarding the integration constants can be partially broken once a complete post- 
Newtonian parameterization is developed in presence of matter. In such a case, the integration constants remain 
constrained by the Boltzmann-Vlasov equation which describes the conservation of matter at these scales [35[ • 

Up to now, the discussion has been developed without any gauge choice. In order to overcome the difficulties related 
to the nonlinearities of calculations, we can work considering some gauge choice obtaining less general solutions for the 
metric entries. A natural choice is represented by the conditions ([20]) which coincide with the standard post-Newtonian 



** In this case the symbol 0[2] is referred to higher than first order contributions the dimensionless coordinate f. 
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gauge 



hjk, ' 



1 



h,j = 0(4), 



h Qk «--h k k , = 0(5) 



(45) 



where accounts for deviations from the Minkowski metric (<7 M „ = rj^ u + h^). In this case the Ricci curvature 
tensor becomes 1 



R 



tt\ h g 



R 



(2) 



R 



rr\h g 



R 



R 



(4) 



(2) 



(46) 



where 



R 



(2) _ rg£] rT +2g£] T 



R 



tt\ hg 2r 

(4) _ r fl ff r r + 2g ff„ +r[gg g <?> rr - g ff ft - g &> r . 



2r 



i? 



(2) _ rg^ rr +2g^ r 



2r 



R {2) = ^(2) = Q 



while the Ricci scalar expressions at the 0(2) and 0(4) orders read 



R 



(2) _ rg™„+2 g ™ r -rg™„-2g™ r 



\hg 



2r 



(47) 



R 



(4) = ^+ggg>jgg^ 
|ftg 2r 



(48) 



The gauge choice does not affect the Christoffel. Thus, by solving the system with the simplification induced by 
the gauge, one obtains 



g (t,r) = l + ^ 



ki , k a log r 



(49) 



where the constants ki, k± are relative to the 0(2) order of approximation, while ki and k^ are related to the 0(4) 
order. The Ricci scalar is zero both at 0(2) and at 0(4) approximation orders. 

Eqs. (|49|) suggest some interesting remarks. It is easy to check that the GR prescriptions are immediately recovered 
for k\ = ki and k2 — k% — 0. The g rr component displays only the second order term, as required by a GR-like 
behavior, while the gtt component shows also the fourth order corrections which determine the second post-Newtonian 
parameter (3 [25| . It has to be stressed here that a full post-Newtonian formalism requires to take into account matter 
in the system (|34j) : the presence of matter links the second and fourth order contributions in the metric coefficients 

m. 



1 We have indicated with the subscript hq the harmonic gauge variables 
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V. THE POST-MINKOWSKIAN APPROXIMATION 



In the previous section we have developed a general analytic procedure to deduce the Newtonian and the post- 
Newtonian limit of /(-R)-gravity in absence of matter or far from matter sources. Here we want to discuss a different 
limit of these theories, pursued when the small velocity assumption is relaxed and only the weak field approximation 
is retained. This situation is related to the Minkowski limit of the underlying gravity theory as well as the discussion 
of the previous section was related to the Newtonian one. In order to develop such an analysis, we can reasonably 
resort to the metric (j4j, considering the gravitational potentials A(t, r) e B(t, r) in the suitable form 



A(t,r) = l + a(t,r) 
B(t,r) = l + 6(t,r) 



(50) 



with a(t,r),b(t,r) <C 1. Let us now perturb the field equations ([5|) considering, again, the Taylor expansion (|33|) for 
a generic f(R) theory. For the vacuum case (T^ v — 0), at the first order with respect to a e b, it is 



' fa = 

/i{4i 3 -|^ (1) }+^ = o 



(51) 



where 



W« - -/ate - r (0) >« - g® gWj® + g®"H$ + In y/=$B® 



(52) 



In this approximation, the Ricci scalar turns out to be zero while the derivatives, in the previous relations, are 
calculated at R = 0. 

Let us now consider the limit for large r, i.e. we study the problem far from the source of the gravitational field. 
In such a case Eqs. (|5ip become 



d'ajt r) _ d 2 b{tr) 
d? 1 dP 







a(t, r) - b(t, r) 



d 2 b(t,r) d 2 a(t.r) _ n d 2 b(t,r) 

9r 2 T g t 2 * Q t 2 



*(*) 



(53) 



where ^(t) is a generic time-dependent function. Eqs. (|53[) are two coupled wave equations in term of the two functions 
a(t, r) and b(t, r). Therefore, we can ask for a wave-like solutions for the gravitational potentials ait, r) and b(t, r) 



a(t, r) = J ^a(uj, k)e t{ut ~ k ^ 
b(t, r) = j <^b(u, fc)e i (" t - fer ) 



(54) 



and substituting these into the (|53p . In order to simplify the calculations, we can set ^(t) = since, as said, this is 
an arbitrary time function. Eqs. (|53|) are satisfied if 



a(uj, k) = b(u>, k) , uu = ±k 

< 

d(tu, k) — 

fi 

where, as before, £ = - L — . In particular, for fi — or fa = one obtains solutions with a dispersion relation u = ±fe. 
6/2 

In other words, for /j ^ (i = 1,2), that is in the case of non-linear f(R), the above dispersion relation suggests 
that massive modes arc in order. In particular, for £ < 0, the mass of the graviton is m grav = — — and, coherently, 
it is obtained for a modified real gravitational potential. As matter of fact, a gravitational potential deviating from 



i_ JL 

1 4fe 2 



b{u,k), uj^±\/k 2 -^ 
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the Newtonian regime induces a massive degree of freedom into the particle spectrum of the gravity sector with 
interesting perspective for the detection and the production of gravitational waves [36|. It has to be remarked that 
the presence of massive gravitons in the wave spectrum of higher order gravity is a well known result since the paper 
of [111 ] . Nevertheless it is our opinion that this issue has been always considered under a negative perspective and has 
been not sufficiently investigated. Furthermore, if £ > 0, even the solution 



a(t, f) = (ao + air)e ± ~ t 



b(t, f) = (bo + bit) cos 



V3f 
2 ' 



O'o + sin 



2 ' 



fig + b'{t 



(56) 



with ao, ai, bo, b\, b' , bi, 6 ', b" constants is admitted. The variables f and t are expressed in units of £ -1 / 2 . In the 
post-Minkowskian approximation, as expected, the gravitational field propagates by means of wave- like solutions. This 
result suggests that investigating the gravitational waves behavior of fourth order gravity can represent an interesting 
issue where a new phenomenology (massive gravitons) has to be seriously taken into account. Besides, such massive 
degrees of freedom could be a realistic and testable candidate for cold dark matter, as discussed in [37j |. 



VI. CONCLUSIONS 



In this paper, we have developed a general analytic approach to deal with the weak field and small velocity limit 
(Newtonian limit) of a generic f(R) gravity theory. The scheme can be adopted also to correctly calculate the post- 
Newtonian parameters of such theories without any redefinition of the degrees of freedom by some scalar field leading 
to the so called O'Hanlon Lagrangian [38|. In fact, considering this latter approach, we get a Brans-Dicke like theory 
with vanishing kinetic term and then the post-Newtonian parameter 7 results 7=1/2 and not 7 ~ 1 as observed. 
This result is misleading in the weak field limit. In the approach presented here, we do not need any change from the 
Jordan to the Einstein frame [3(1 H^|. Apart the possible shortcomings related to non-correct changes of variables, 
any f(R) theory can be rewritten as a scalar-tensor one or an ideal fluid, as shown in [4(| [42j ■ In those papers, it 
has been demonstrated that such different representations give rise to physically non-equivalent theories and then also 
the Newtonian and post-Newtonian approximations have to be handled very carefully because the results could not 
be equivalent. In fact, the further geometric degrees of freedom of f(R) gravity (with respect to GR), the scalar field 
and the ideal fluid have weak field behaviors strictly depending on the adopted gauge which could not be equivalent 
or difficult to compare. In order to circumvent these possible sources of shortcomings, one should states the frame 
(Jordan or Einstein) at the very beginning and then remain in such a frame along all the calculations up to the final 
results. Adopting this procedure, arbitrary limits and non-compatible results should be avoided. 

In this paper, we have considered the Taylor expansion of a generic f(R) theory, obtaining general solutions in 
term of the metric coefficients up to the third order of approximation when matter is neglected. In particular, the 
solution relative to the g u metric component gives the gravitational potential which is corrected with respect to 
the Newtonian one of f(R) = R- The general gravitational potential is given by a couple of Yukawa-like terms, 
combined with the Newtonian potential, which is effectively achieved at small distances. In relation to the sign of the 
characteristic coefficients entering the g u component, one can obtain real or complex solutions. In both cases, the 
resulting gravitational potential has physical meanings. This degeneracy could be removed once standard matter is 
introduced into dynamics. 

The complete analysis allows to obtain direct information on the post-Newtonian formalism: the post-Newtonian 
parameters can be fully characterized considering the integration constants in the gravitational potential. Nevertheless 
this study is beyond the aim of this paper and will be developed in a forthcoming research project. 

Furthermore, it has been shown that the Birkhoff theorem is not a general result for /(i?)-gravity. This is a 
fundamental difference between GR and fourth order gravity. While in GR a spherically symmetric solution is, in any 
case, stationary and static, here time-dependent evolution can be achieved depending on the order of perturbations. 

Finally, we have discussed the differences between the post-Newtonian and the post-Minkoskian limit in f(R) 
gravity. The main result of such an investigation is the presence of massive degrees of freedom in the spectrum of 
gravitational waves which are strictly related to the modifications occurring into the gravitational potential. This 
occurrence could constitute an interesting opportunity for the detection and investigation of gravitational waves. 



15 



VII. APPENDIX 

In this appendix, we show that the harmonic gauge can be suitably reduced to the form (|20|) . Such a gauge is 
usually characterized by the condition g aT T^ T = 0. For /i = one has 

2 9 -Tl T ^ 9 K ~^T+9 {2 C° = 0, (57) 

and for fi = i 

2grKr « 5 (2) o + 2ff (2) , m - 9 {2) m = . (58) 
Differentiating Eq. ([57]) with respect to a; , x 3 and ([5"5]) and with respect to one obtains 

9 0,00 Z 5 0,0m ' 9 m.OO — U > l 0y J 

(2)° _ 2o ( 3 ) m W 2 ) m -o rem 

9^ + 29^-9^ = 0. (61) 
On the other side, combining Eq.(|60p and Eq. ([6l"]) , we get 

(2) m _ fl (2) m „(3) m n f g 2 % 

9 m,0i 9 i,0m f 0,mi — u • l uz 7 

Finally, differentiating Eq. (|5"8|) with respect to a;- 7 , one has : 

9^+29^-9^=0 (63) 
and redefining indexes as j — > i, i — > j since these are mute indexes, we get 

9^+2^-9^=0- (64) 
Combining Eq. ([63]) and Eq. ([64|) . we obtain 

(2)° +0 (2) m , (2) m _ (2) m _ f651 
The relations ([59]). ([62]) , ([65]) guarantee the viability of ([20]) . 
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